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Summary 

The  infrared  region  of  the  electromagnetic  spectrum,  which  covers  wavelengths  from  1  to  100  /./m 
has  many  applications  in  chemical  and  biological  sensing,  thermal  and  night  imaging  etc.  To  fulfill 
these  applications  it  is  necessary  to  have  sources  and  detectors  of  infrared  radiation.  The  latter  are 
usually  semiconductor  or  solid  state  devices  which  have  a  narrow  or  broadband  response  at  a  fixed 
wavelength.  Thus  a  device  often  has  to  be  designed  and  constructed  for  a  particularly  wavelength 
to  detect  a  particular  compound.  This  project  has  been  a  feasibility  study  into  the  use  of  magnetic 
semiconductors  to  make  a  tunable  infrared  detector.  Such  a  device  could  be  tuned  with  a  magnetic 
field  to  different  absorption  lines  of  different  molecules  and  hence  provide  spectral  information  for 
easier  species  identification.  It  could  also  be  used  for  hyperspectral  imaging  to  provide  greater 
resolution  or  adaptable  night  vision/thermal  imaging  to  suit  the  particular  locality  or  environmental 
conditions.  The  general  principle  of  using  magnetic  semiconductors  for  tunable  quantum  well  energy 
structures  could  be  extended  to  tunable  infrared  sources.  As  far  as  the  author  is  aware  no  such  solid 
state  tunable  infrared  technology  is  available. 

The  technical  aspects  of  the  project  are  described  in  the  attached  preprint  which  has  just  been 
submitted  to  the  Journal  of  Applied  Physics  which  is  published  by  the  American  Institute  of  Physics. 


Roadmap  for  realisation  of  tunable  infrared  detectors 

The  central  aspect  of  the  technology  is  the  growth  of  the  magnetic  semiconductor  structures.  Leeds 
itself  is  a  semiconductor  growth  centre  for  GaAs-based  materials  and  as  such  we  have  excellent 
relationships  with  other  growth  centres.  There  are  several  groups  around  the  world  who  could,  given 
the  exact  layer  specification,  grow  this  material.  The  remainder  of  the  technology  could  be  developed 
in-house  in  Leeds.  This  would  include  the  initial  semiconductor  wafer  design,  the  device  design,  the 
device  fabrication  and  the  magneto-infrared  characterisation. 
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Abstract 

A  theoretical  model  and  a  design  of  a  magnetic  field  tunable  CdMnTe/CdMgTe  terahertz  quan¬ 
tum  well  infrared  photodetector  are  presented.  The  energy  levels  and  the  corresponding  wave- 
functions  were  computed  from  the  envelope  function  Schrodinger  equation  using  the  effective  mass 
approximation  and  accounting  for  Landau  quantization  and  the  giant  Zeeman  effect  induced  by 
magnetic  confinement.  The  electron  dynamics  were  modeled  within  the  self-consistent  coupled  rate 
equations  approach,  with  all  relevant  electron-longitudinal  optical  phonon  and  electron-longitudinal 
acoustic  phonon  scattering  included.  A  perpendicular  magnetic  field  varying  between  0  T  and  5  T, 
at  a  temperature  of  1.5  K,  was  found  to  enable  a  large  shift  of  the  detection  energy,  yielding  a 
tuning  range  between  24.1  meV  and  34.3  meV,  equivalent  to  51.4  gm  to  36.1  gm  wavelengths.  For 
magnetic  fields  between  1  T  and  5  T,  when  the  electron  population  of  the  QWIP  is  spin-polarized, 
a  reasonably  low  dark  current  of  <  1.4  x  10~2  A/cm2  and  a  large  responsivity  of  0.36  —  0.64  A/W 
are  predicted. 
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I.  INTRODUCTION 


In  recent  years  there  has  been  an  increasing  interest,  and  a  considerable  experimental 
and  theoretical  activity  focused  on  dilute  magnetic  semiconductors  (DMSs)  whose  uncon¬ 
ventional  properties  offer  a  great  prospect  for  developing  a  wide  range  of  novel  devices1-6. 
The  key  reason  for  employing  them  in  both  basic  and  application-oriented  research  is  the 
control  of  spin  splitting  that  can  be  achieved  by  varying  the  material  composition,  an  ex¬ 
ternal  magnetic  field,  temperature  and/or  via  quantum  confinement7,8.  These  extensive 
studies  have  been  conducted  mostly  in  the  area  of  spintronics,  towards  realization  of  spin- 
based  devices1-4.  Some  of  the  research  activity  has  been  oriented  into  the  investigation  of 
the  spin  splitting  effects  on  interband  transitions  and  the  applications  of  DMSs  in  interband 
devices9-12.  On  the  contrary,  very  few  studies  aiming  towards  novel  intersubband  devices 
based  on  DMSs  have  been  reported13. 

In  spite  of  significantly  improved  performance  and  technological  advances  of  quantum 
well  infrared  photodetectors  (QWIPs)14  and  quantum  cascade  lasers  (QCLs)15,  they  are  still 
attracting  a  considerable  amount  of  research  attention.  Different  designs  of  QWIPs16-19,  the¬ 
oretical  models  of  electron  transport20-23  and  various  applications24,25  have  been  proposed. 
QCL  operation  has  been  recently  demonstrated  in  the  far-infrared  range  up  to  160/jrn26  29 . 
The  development  of  terahertz  (THz)  QCLs  has  motivated  the  research  leading  to  devel¬ 
opment  of  THz  QWIPs30-32.  Recently,  the  experimental  observations  of  greatly  improved 
GaAs/AlGaAs  QCL  performance  under  an  applied  magnetic  field,  in  terms  of  threshold 
current  and  its  temperature  dependence,  have  been  reported29,33,34. 

However,  one  of  the  much  sought  after  improvements  of  THz  QWIPs  comprising  doped 
symmetric  rectangular  quantum  wells  is  external  tunability.  Magnetic  field  appears  to  offer 
great  potential  for  tuning,  if  applied  on  structures  based  on  DMSs.  Due  to  the  manipulation 
of  the  electronic  structure  by  varying  magnetic  field,  one  can  change  the  transition  energies 
of  the  absorbed  radiation13.  Variations  of  a  few  meV  which  are  achievable  in  CdMnTe  (a 
DMS  compound)  imply  a  large  tunability  on  a  relative  scale  if  the  transition  energy  is  in 
the  THz  range. 

The  aim  of  the  present  work  is  to  explore  the  design  of  a  magnetically  tunable  THz  QWIP 
comprising  CdMnTe/CdMgTe  materials.  Additionally,  introducing  a  new  material  system 
for  QWIPs  offers  a  possibility  to  access  the  GaAs  Reststrahlen  region  of  34  —  36  meV.  A 
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model  of  electron  transport  in  QWIPs  in  a  magnetic  field  is  then  required  in  order  to  predict 
the  QWIP  output  characteristics. 

The  organization  of  the  paper  is  as  follows.  In  Sec.  II  a  theoretical  model  of  QWIPs 
in  a  magnetic  field  is  presented.  The  calculation  of  the  electronic  structure  (energy  levels 
and  wavefunctions)  of  DMS  quantum  structures  is  described  in  Sec.  IIA.  The  interaction 
of  electrons  with  electromagnetic  radiation  in  an  external  magnetic  held  in  the  dipole  ap¬ 
proximation  is  reviewed  in  Sec.  IIB.  The  scattering  rates  between  Landau  levels  due  to  the 
interaction  with  optical  and  acoustic  phonons  are  given  in  Sec.  IIC.  Eventually,  the  rate 
equation  model  of  QWIPs  subjected  to  a  magnetic  held  is  presented  in  Sec.  IID.  The  details 
of  design  and  the  results  are  presented  in  Sec.  III.  The  tuning  range  of  the  QWIP  in  a  mag¬ 
netic  held  up  to  5  T  and  the  variation  of  the  electron  distribution  are  presented  in  Sec.  IIIA 
and  Sec.  IIIB,  respectively.  The  dark  current  and  responsivity  versus  magnetic  and  electric 
held  dependences  are  discussed  in  Sec.  IIIC  and  Sec.  HID. 


II.  THEORETICAL  CONSIDERATIONS 

A.  Electronic  structure  of  DMS  quantum  wells 

The  unique  properties  of  DMSs  arise  from  their  band  structure  with  two  distinct  elec¬ 
tronic  subsystems:  delocalized,  band  electrons  and  magnetic  impurity  electrons  with  mag¬ 
netic  moments  localized  in  ionic  open  3 d  (or  4/)  shell.  The  mobile  electrons  determine 
the  electrical  and  optical  properties  of  DMSs,  while  the  localized  magnetic  moments  are 
responsible  for  their  magnetic  properties.  However,  the  unique  magneto-optical  properties 
of  DMSs  result  from  the  strong  spin-dependent  sp-d(f)  exchange  interactions  between  these 
two  subsystems5”8. 

Application  of  a  perpendicular  magnetic  held  to  nonmagnetic  semiconductor  quantum 
wells  splits  the  in-plane  continuum  of  quantized  subbands  into  Landau  levels,  each  subband 
producing  a  set  of  Landau  levels,  described  by  a  Landau  index  and  additionally  by  spin 
index  (gyromagnetic  spin  splitting  within  Landau  levels)35.  Both  the  Landau  level  energies 
and  their  separation  depend  linearly  on  the  magnetic  held.  If  the  structure  includes  a  DMS, 
the  conduction  band  edge  varies  with  the  magnetic  held  locally  due  to  the  sp-d  interaction, 
and  hence  modifies  the  potential  profile.  The  shift  of  the  potential  (in  the  magnetic  layers 


3 


and  near  the  interfaces)  is  opposite  for  the  two  spins,  and  consequently  they  experience  two 
different  potential  profiles,  which  are  both  a  function  of  the  magnetic  held  5'9>10>36.  This  spin 
dependent  variation  of  the  confinement  energy  is  particularly  significant  in  the  vicinity  of 
the  interfaces10,36  and  in  turn  leads  to  different,  and  field-dependent,  spectra  for  the  two 
spins.  This  phenomenon  is  known  as  the  giant  Zeeman  effect.  Since  the  spin  and  Landau 
level  index  are  conserved  in  optical  transitions  caused  by  ^-polarized  light,  the  transition 
energies  will  vary  with  the  magnetic  held.  This  translates  into  tunability  of  the  intersubband 
transition  energies  by  varying  the  magnetic  held. 

The  total  Hamiltonian  of  an  electron  in  DMSs,  in  addition  to  the  conventional  Hamil¬ 
tonian  for  nonmagnetic  semiconductor  structures,  contains  the  spin-dependent  potential 
induced  by  the  magnetic  held,  given  by  an  empirical  expression5,9,10,36 

UDms  =  ±^aN0x  (Sz) ,  (1) 


where  “+”  holds  for  spin-up  and  ”  for  spin-down  electrons,  a  is  the  sp-d  exchange  integral 
for  the  conduction  band,  Nq  is  the  density  of  cations,  (Sz)  is  the  thermal  average  of  the  Mn2+ 
spin,  which,  for  paramagnetic  materials,  amounts  to 


(S,) 


5  5  g*pBB 

2  5/2  [2k(T  +  T0) 


(2) 


where  g*  is  the  Lande  factor  of  the  conduction  band  electrons  (here  assumed  constant 
throughout  the  structure,  g*(z)  ~  2),  /iB  is  the  Bohr  magneton,  B  is  the  magnetic  held 
along  the  z  axis  (the  direction  of  confinement),  k  is  the  Boltzmann  constant,  T  the  temper¬ 
ature  and  H5/2  is  the  Brillouin  function 

w)  —  H  coth  ^  coth  .  (3) 


Phenomenological  fitting  parameters  x  and  T0  represent  the  effective  manganese  concentra¬ 
tion  and  the  effective  temperature5,9,10,36  and  both  depend  on  the  manganese  concentration 
and  the  position  on  the  z  axis  if  the  influence  of  a  graded  nonDMS/DMS  interface  and 
the  enhanced  magnetization  in  its  vicinity  is  accounted  for.  A  theoretical  method  mapping 
the  detailed  interface  profile  was  used  to  obtain  spatial  and  manganese  concentration  de¬ 
pendences  of  x  and  T036.  The  band  offset  was  calculated  accounting  for  the  variation  of 
composition  profile  of  Mn  ions  in  the  magnetic  barriers  and  at  the  interface. 
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The  electronic  structure  of  DMS  quantum  wells  was  found  by  solving  Schrodinger  equa¬ 
tion  for  envelope  functions  within  the  effective  mass  approximation 

h2  d  (  1  d^(z)\  ,  r^._,  ,  i_  ,  1  =  V1  ,,, 


2  dz  \m*(z)  dz 


+  Ec(z)  +  eFzT  ^9*(z)l^BB  ± -aN0x  (Sz)  ip(z) 


=  \E 


J  +  2 


1\  heB 


where  h  is  Planck  constant,  e  is  the  electron  charge,  rn*  is  the  electron  effective  mass.  The 
second  term  Ec(z)  on  the  left  side  of  Eq.  4  is  the  band  offset,  the  third  term  is  the  potential 
due  to  an  applied  electric  field  F,  the  fourth  is  the  potential  stemming  from  the  gyromagnetic 
spin  splitting  of  Landau  levels  (“— ”  for  spin- up  and  “+”  spin-down  electrons),  and  the  fifth 
is  the  spin-dependent  potential  of  DMS.  In  further  considerations  we  will  use  a  shorthand 
subscript  i  to  denote  the  jtth  Landau  level  (LL)  of  the  m,th  state  (subband)  with  spin  s*, 
i.e.  i  =  | rrii,ji,Si).  Its  energy  is 


Ei  —  E\mi  —  En ^ ' ji  +2^ 


where  E ^  is  the  energy  of  the  state  rnt  with  spin  Sj,  and  ujc  =  eB/m*  is  the  cyclotron 
frequency.  The  wavefunctions  of  all  the  LLs  of  the  same  spin,  associated  to  the  same  state, 
have  the  same  form  'ifti(z)  =  'ijj^.(z). 


B.  Interaction  of  electrons  with  electromagnetic  field 


In  optical  transitions  the  Landau  index  and  spin  are  conserved,  and  absorption  occurs 
only  on  transitions  between  LLs  associated  with  different  states.  The  fractional  absorption 
of  such  a  transition  is13,3' 

Ai,f  =  ^M^mfS{Ef  -  ^  -  hu)(nf  -  m),  (6) 

lltQC  J 

where  i  =  \rrii,j,s)  and  /  =  \rrif,j,s)  are  LLs  stemming  from  different  states  (m,  ^  m/ ), 
Ei  and  Ef  are  their  energies  with  the  corresponding  wavefunctions  'ipm,  anfl  = 

I '^Sm*f{z)z'll;mi{z)  dz  is  fhe  dipole  matrix  element,  nt  and  rif  are  the  electron  sheet  densities, 
n  is  the  refraction  index,  eo  is  the  permittivity  of  vacuum,  c  is  the  speed  of  light  in  vacuum, 


and  uj  is  the  frequency  of  incident  radiation.  The  fractional  absorption  may  be  written  as 
Aij  =  CTijirii  —  rtf)  with  the  cross-section 


alJ  =  —M^  6(Ef-El-hu;), 

TICqC  j 
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independent  of  the  electron  concentration.  The  total  absorption  between  states  rril  and  m/ 
of  the  same  spin  s  is  a  sum  of  contributions  of  all  LLs  stemming  from  them 


7re2ca 

neve 


Mf  2 


mi,rrif 


S(K,  -  E, 


mi 


As 


X)  Ahf 

hf,3i=3f=3,Si=Sj=s 


huj)  X  (nf~ni)i 

i,f,3i=3f=3,Si=Sj=s 


(8) 


where  Ejh  —  E^  —  Ei  —  Ef.  State  broadening  was  modeled  by  a  Lorentzian  distribution 
with  the  transition  linewidth  (FWHM)  T. 


C.  Interaction  of  electrons  with  phonons 


Within  the  frame  of  a  bulk  phonon  model,  the  electron-longitudinal  optical  (LO)  phonon 
scattering  rate  between  LLs  i  =  rnt ,  jt ,  s)  and  /  =  \rrif,jf,s)  if  jx  <  jf  reads38 


WLO{i,f)  = 


e2ijJLO 

2'K€ri 


noi^Lo)  +  2^2 


x 


Qxy 


Hji,jf(Qxy )  dqxy 


5(Ei  -  Ef  ±  hujLo ) 

2 

-  dqz, 


r*s 

^  mi,rrif 


( qz 


qly  +  ql 


(9) 


where  Ulo  is  optical  phonon  frequency,  qxy  (qxy  =  qxy  | )  is  the  ^-component  and  qz  is 
the  ^-component  of  the  phonon  wave  vector  q  =  (qxy,  qz%  the  constant  ep  is  dehned  as 
1  =  e"1  —  e^1  where  and  es  are  high-frequency  and  static  permittivity,  respectively, 
no {ojlo)  is  the  Bose-Einstein  factor.  To  calculate  the  scattering  rate  if  jt  >  jf,  one  swaps 
these  indices.  The  upper  sign  in  the  (n0(uoLO)  +  \  T  |)  term  and  in  the  energy  conservation 
law  (the  6  function)  holds  for  absorption,  and  the  lower  sign  for  emission. 

The  form-factor  Gs  m  (qz )  is  dehned  as 

1 1  vj 


G 


rrii,mf 


(?«)  =  /  dz. 


(10) 


The  lateral  overlap  integral  Hju3f(qxy )  is  given  in  the  analytical  form 


H3iJf(Qxy)  I  =e  9 


'I  /  2  \3f-3i 

2  -  p-qly/W^hL  (  Si E  1 


jf'-  V2/52 


j3f  Ji 
L3f 


y.xy 

2^ 


(11) 


m*ojc  _ 

h 


where  f5  =  \J 


and  L*j(x)  is  the  associated  Laguerre  polynomial.  The  delta 
function  was  replaced  by  a  Gaussian39  with  the  transition  linewidth  (FWHM)  equal  to  T  if 
rrii  S1  mj  and  j%  =  j3 ,  or  equal  to  a0\^B\/2  In  2  if  m,;  =  m3  and  jt  S1  jj,  or  1°  \f F2  +  2  ln2oQ.E> 
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if  mi  7^  nij  and  jl  ^  jj.  The  parameter  ctq  is  inversely  proportional  to  the  effective  mass,  so 
from  its  value  for  GaAs39  (1  meV/T1/2),  one  may  estimate  the  value  for  CdMnTe. 

The  transition  of  an  electron  from  the  ith  LL  ( i  =  \ rrii,ji,s))  with  larger  energy  than 
that  of  the  final  / th  LL  (/  =  \rrif,jf,s))  via  an  interaction  with  longitudinal  acoustic  (LA) 
phonons  in  a  magnetic  field  is  given  by40 


WAC(i,f) 


D\(E,  -  E,f 


Ei~Ef 

e  kT 


npVgfi'5 


-  1 


fQzn 


X 


S~ IS 


( Q.Z )  Eji,jf(Qxyo)  d-5; 
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where 


QxyO 


\ 


Ei  ~  Ef 
Ttv „ 


(13) 


Here  Da  is  the  deformation  potential,  p  is  the  density  of  the  material,  vs  the  sound  velocity, 
Gmi,mf(qZ )  anfl  Hjujf(qxy )  are  the  form  factor  and  the  lateral  overlap  integral  for  clectron-LO 
phonon  scattering,  and 


Qzmax 


Ei  ~  Ef 
hvs 


(14) 


If  Ei  <  Ef,  the  electro n-L A  scattering  is  calculated  from  a  similar  expression 

D\(Et  -  Ef)2  1 


WAC(i,f)  = 


TTpVgk6 


e  kT  —  1 


f  Qzmax 

X  7o  \E'mi,mf(qz)\  {QxyO)  \  dqz-  (15) 

Generally,  electron-LA  scattering  is  less  significant  than  clectron-LO  scattering. 

Landau  quantization  and  the  giant  Zeeman  effect  have  a  significant  impact  on  scatter¬ 
ing  processes,  enhancing  or  inhibiting  them  depending  on  the  inter-LL  energy  separation. 
The  clectron-LO  phonon  interaction  between  a  ground  state  LL  and  a  continuum  state  LL 
increases  considerably  as  their  energy  difference  approaches  one  LO  phonon  energy  (Elq)- 
The  scattering  reaches  a  maximum  when  the  difference  becomes  equal  to  Elq  and  then  the 
corresponding  LLs  are  in  resonance.  A  similar  conclusion  applies  for  electron-LA  phonon 
scattering,  with  the  distinction  that  resonance  occurs  when  two  LLs  have  very  similar  ener¬ 
gies. 
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D.  Modeling  of  QWIPs  in  a  magnetic  field 


A  quantum  mechanical  model  of  electron  dynamics  in  QWIPs  in  a  magnetic  field,  ac¬ 
counting  for  Landau  quantization  and  the  giant  Zeeman  effect,  was  developed  to  describe 
the  electron  transport  processes.  Consider  a  QWIP  with  a  large  number  of  periods  in  an 
externally  applied  electric  and  magnetic  field.  The  population  dynamics  of  LLs  can  be 
obtained  from  the  system  of  nonlinear  rate  equations  in  the  steady  state21,41 
dn  nc  nc  nc 

-*r  =  0  =  53  riiWij—nf  WU  +  C{$,ni,nf),  (16) 

al  i=U^f  »=lMf 

where  i  and  /  run  over  all  Nc  LLs  in  the  cascade,  in  all  of  its  periods,  nt  is  the  electron 

concentration  of  the  ith  LL  and  IF,;,/  is  the  total  scattering  rate  from  LL  i  into  LL  / 

W,J  =  W.jll  -  fMEf)],  (17) 


where  W ij  is  the  scattering  rate  from  LL  i  into  LL  /,  independent  of  the  electron  distribu¬ 
tion,  and  l  —  fFo(Ef)  accounts  for  Pauli  exclusion,  with  fpD  being  the  Fermi  Dirac  function. 
The  electron  distribution  over  all  LLs  associated  with  some  state  is  assumed  to  be  Fermi 
Dirac-like,  with  electron  temperature  equal  to  the  lattice  temperature,  but  with  different 
quasi-Fermi  level  EFm.  for  each  state  m*,  i.e. 


IfD  (.Ei)  Bi-Bp, m.  > 

e  kT  —  1 

The  electron  concentration  on  the  ith  LL,  with  energy  Ell  is 


(18) 


ni  —  n\ m,i,ji,Si) 


eB 
2n  h 


fpD^Ei). 


(19) 


The  last  sum  in  Eq.  16  is  the  contribution  of  intersubband  absorption  to  electron  transition 
rates  with21 

C(<F,  m,  nf)  =  Aif(ni ,,  nf)$,  (20) 

where  Atf  is  the  fractional  absorption  on  the  i  — >  f  transition,  and  <F  the  optical  flux. 

Electron-LO  phonon  and  electron-LA  phonon  scattering  are  included  in  the  model  as 
the  main  scattering  mechanisms.  Ionized  impurity  scattering  is  considered  negligible  as  the 
doping  density  is  assumed  low.  Electron-electron  scattering  between  the  ground  state  LLs, 
populated  by  the  majority  of  carriers,  and  the  LLs  in  the  continuum  is  neglected  because  of 


the  relatively  large  energy  separation,  as  well  as  between  any  two  of  the  continuum  state  LLs 
because  of  their  low  population.  Spin-flip  scattering  is  in  general  less  significant  than  phonon 
scattering42  and  was  neglected.  Since  the  scattering  between  electrons  of  different  spin  is 
not  included  in  the  model,  the  rate  equation  system,  which  breaks  into  two  independent 
subsystems,  one  for  either  spin  orientation,  together  with  the  particle  conservation  law, 
is  under  determined.  Therefore,  one  needs  to  calculate  the  total  concentrations  of  spin-up 
and  spin-down  electrons.  Since  the  perturbation  of  the  system  induced  by  an  external 
electric  held  is  relatively  small,  a  quasiequilibrium  state  of  the  system,  with  a  unique  quasi- 
Fermi  energy  Ep,  is  assumed  (only  for  calculating  the  total  spin-up  and  spin-down  electron 
concentrations).  Inserting  the  quasi- Fermi  energy  Ep  into  Eq.  18,  replacing  Eq.  18  in  Eq.  19 
and  summing  over  all  the  LLs  one  may  obtain  Ep,  then  End  the  concentrations  rq  from 
Eq.  19  and  sum  those  of  the  same  spin.  Hence,  each  electron  spin  subsystem  is  effectively 
described  by  an  independent  rate  equation  system  and  a  particle  conservation  law.  In  this 
paper  the  simulations  were  carried  out  in  the  magnetic  held  range  between  1  T  and  5  T, 
where  the  population  of  spin- up  states  was  found  to  be  negligible  (given  in  Sec.  IIIB).  Hence, 
only  spin-down  electron  transitions  were  taken  into  account  in  the  analysis,  using  only  the 
rate  equations  and  the  particle  conservation  law  for  the  spin-down  electrons. 

Due  to  the  wavefunction  localization  properties,  each  state  and  the  LLs  originating  from 
it  can  be  associated  to  one  of  the  periods  of  the  QWIP.  We  assume  a  globally  linear  variation 
of  the  conduction  band  potential  and  apply  the  quasiperiodicity  principle,  i.e.  if  ^(z)  is  the 
wavefunction  of  a  LL  of  energy  Et,  then  -0j( z  —  D )  is  the  wavefunction  of  a  LL  of  energy  E)  — 
AH,  where  D  is  a  period  length  and  AH  is  the  potential  drop  across  a  period.  Therefore,  by 
determining  all  the  LLs  (their  energies  and  wavefunctions)  belonging  to  one  period  (referred 
to  as  the  central),  one  may  obtain  all  the  LLs  in  the  whole  structure.  It  is  considered  that 
a  state  and  its  LLs  belong  to  a  period  if  that  state  has  a  better  overlap  with  the  ground 
state  of  that  period  than  with  the  ground  states  of  the  neighboring  periods  (overlap  is  here 
defined  as  the  integral  of  their  squared  wavefunctions  along  the  structure).  An  identical 
set  of  N  states  and  their  LLs  with  potentially  considerable  population  are  thus  attributed 
to  each  period  (high  energy  states  are  practically  unpopulated).  This  discretization  of  the 
continuum  spectrum  gives  an  accurate  description  of  the  continuum  if  addition  of  new  states 
and  their  LLs  from  the  same  basis  does  not  affect  the  values  of  physical  observables  i.e. 
current.  A  LL  assigned  to  the  central  period  is  labeled  as  the  ith,  ( i  =  1,  ...IV),  and  a 
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LL  assigned  to  its  kt\i  nearest  neighbor  as  the  i  +  kN th,  (i  =  1  where  k  >  0  for 

right  and  k  <  0  for  left  neighbors.  Since  the  interaction  between  LLs  belonging  to  distant 
periods  is  negligible  clue  to  the  small  wavefunction  overlap,  the  ‘tight-binding’  description 
is  introduced,  by  including  the  interaction  with  only  P  nearest  neighbors  on  both  sides. 
Assuming  an  identical  electron  distribution  in  each  period  (n,  =  ni+kN),  one  may  account 
explicitly  for  distributions  in  the  central  period,  and  evaluate  the  scattering  rates  within  it 
and  between  any  two  LLs  belonging  to  the  central  period  and  any  of  its  P  nearest  neighbors 
on  either  side.  Furthermore,  the  scattering  rates  are  shift  invariant  ( =  Wi+kN,j+kN, 
Wi-kN,j  =  t'ldy+fcv),  therefore  reducing  the  calculation  of  the  scattering  rates  by  taking  into 
account  only  the  central  period  and  its  P  right  neighbors.  Introducing  olb  =  A h  in  Eq.  19 
and  from  Aij  =  —  rif )  and  Eq.  7,  we  get  a  system  of  nonlinear  equations 

(jrif  n  _  N 

—jr  =  0  =  X!  ni\yVi,fO-  -  aB'nf )  +  -  nf  E  iW fA1  -  aB'rii )  + 

at  i  l.</f  i-  1  •'// 

P  N 

+  E  E  {rii[(yv i,f+kN  +  Wi+kN,f){  1  —  OlB'rif )  +  (cTjj+jfcjv  +  <Ji+kN,f)&] 

k=l  i=l,i^zf 

—nf[(Wf+kN,i  +  W//,j+fcAf)(l  —  OlB'rii)  +  ((Jf+kN,i  +  C/.j+fcTv)®]}  (21) 


where  i  and  /  run  over  all  LLs  in  a  period.  Within  the  ‘tight-binding’  description,  the 
interaction  with  two  nearest  neighbors  (P  =  2)  is  considered. 

The  current  density  is  calculated  by  taking  into  account  all  electrons  crossing  some  ref¬ 
erence  plane  e.g.  the  interface  between  the  central  period  and  the  adjacent  right  period. 
Effectively,  this  is  done  by  subtracting  the  current  density  component  due  to  electrons  scat¬ 
tering  through  the  reference  plane  in  the  direction  of  the  potential  drop  from  the  component 
due  to  electrons  scattering  in  the  opposite  direction 

P  N  N 

J  =  EEE  krii  [Wij+kN  —  Wi+fciV,/  +  (o'iJ+kN  —  C?i+kN,f)  *£] 

k=  1  i= 1  /=  1 


P  N  N 

=  E  E  E  kni 

k= 1  i=  1  /= 1 


(Wi 


i,f+kN  ~  W i+k N,f)  (1  _  OisUf)  +  {<JiJ+kN  ~  O’i+kNj)  $ 


(22) 


The  responsivity  of  a  QWIP  is  evaluated  from 


J($)  -  J($  =  0) 
{hc/\)$ 


where  A  is  the  detection  wavelength. 


(23) 


10 


III.  NUMERICAL  RESULTS 


The  design  targets  at  simultaneously  achieving  maximum  tunability  and  large  oscillator 
strength  of  transitions  relevant  for  the  device  operation.  These  clearly  are  the  transitions 
of  spin-down  electrons  from  the  ground  state  to  resonant  states  near  the  barrier  top43.  A 
good  design  should  maintain  the  above  features  throughout  the  range  of  magnetic  fields, 
from  zero  to  where  the  sp-d  interaction  saturates. 

The  design  cycle  proceeded  by  varying  the  thickness  and  the  composition  of  the  magnetic 
Cdi_xMnxTe  wells  and  the  nonmagnetic  Cdi_yMgyTe  barriers.  Also,  special  care  was  taken 
to  minimize  the  scattering  rate  from  the  ground  state  to  the  continuum,  in  order  to  minimize 
the  dark  current.  This  was  done  by  tailoring  the  transition  energies  to  be  approximately  one 
linewidth  larger  than  the  LO-phonon  energy  of  the  CdMnTe/CdMgTe  system  (21  meV).  The 
linewidth  was  assumed  to  be  30%  of  the  transition  energy  of  maximal  strength,  as  is  usually 
done  in  bound-continuum  QWIPs44.  Solving  the  Schrodinger  equation  was  implemented  by 
a  finite- difference  method.  Many  body  effects  were  not  taken  into  account  in  the  calculation 
due  to  the  low  doping  density.  The  material  parameters  for  CdMnTe/CdMgTe  in  the  electron 
structure  computation  were  taken  from  Refs. 45-4 '.  The  parameters  used  in  the  scattering 
rate  calculation  were  taken  from  Refs.46,48. 

A.  Tunability 

The  designed  QWIP  consists  of  a  series  of  8.64  nm  Cdo.89M1io.nTe  magnetic  wells  and 
86.4  nm  Cdo.sssMgo.usTe  barriers.  The  doping  density  of  the  wells  was  chosen  to  be 
10ncm-2,  to  satisfy  the  requirement  for  achieving  relatively  high  absorption  in  the  THz 
range  at  low  temperatures31  (~  0.04%).  The  temperature  was  set  to  1.5  K,  when  the  giant 
Zeeman  effect  is  more  pronounced  at  lower  magnetic  fields.  The  conduction  band  edge  po¬ 
tential  profile  of  the  QWIP  in  magnetic  and  electric  fields  for  spin-down  electrons,  together 
with  the  energies  of  the  relevant  states  assigned  to  one  period  is  given  on  the  left-hand  side 
of  Fig.  1.  The  modification  of  the  conduction  band  edge  for  spin-up  and  spin  down-electrons 
induced  by  the  magnetic  field  is  shown  in  the  right-hand  side  of  Fig.  1. 

With  an  applied  electric  field  of  1  kV/cm  and  no  magnetic  field,  the  intersubband  transi¬ 
tion  energy  corresponding  to  the  maximum  oscillator  strength  is  26.9  meV.  It  increases  with 
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magnetic  field  to  34.3  meV  for  spin-down  and  decreases  to  19.5  meV  for  spin-up  transitions 
until  the  sp-d  exchange  interaction  saturates  at  B  —  5  T,  Fig.  2.  With  this  design,  the 
energy  range  of  19.5  —  24.1  meV  is  not  accessible  because  of  the  small  population  of  spin-up 
electrons  for  magnetic  fields  of  1  —  5  T  (explained  below),  effectively  limiting  the  tuning 
range  to  24.1  —  34.3  meV. 

The  contribution  of  other  continuum  states  to  the  responsivity  slightly  shifts  the  peak 
detection  energy  upwards.  For  magnetic  fields  of  1  —  5  T,  for  spin-down  transitions,  the 
peak  detection  energy  differs  from  the  transition  energy  of  maximal  strength  up  to  1.5  meV, 
maintaining  the  tunability  obtained  from  the  consideration  of  the  maximal  strength  tran¬ 
sitions.  ft  is  expected  that  this  conclusion  also  applies  for  both  the  spin-up  and  spin-down 
transitions  in  magnetic  fields  of  0  —  1  T. 

B.  Population  of  Landau  levels 

The  population  of  LLs,  especially  those  stemming  from  the  ground  state,  determines  the 
QWIP  output  characteristics,  such  as  current  and  responsivity.  Therefore,  the  discussion  of 
the  LL  population  is  given  first. 

The  energies  of  several  lowest  index  LLs  originating  from  the  ground  and  the  first  contin¬ 
uum  state  of  both  the  spin-up  and  spin-down  subsystems,  and  the  quasi-Fermi  level  EF  vs. 
magnetic  field  dependences  are  given  in  Fig.  3.  The  inset  of  Fig.  4  shows  the  distribution 
of  spin-up  and  spin-down  electrons  in  a  magnetic  field  from  1  T  to  5  T.  At  magnetic  fields 
above  1  T  the  LL  configuration  is  such  that  all  the  LLs  of  spin-up  electrons  are  above  EF. 
Given  that  the  temperature  is  near  zero,  the  population  of  these  levels  is  extremely  small, 
and  so  is  their  contribution  to  the  absorption  and  scattering  processes.  For  spin-down  elec¬ 
trons  some  low  index  LLs  of  the  ground  state  are  below  EF  and  hence  most  of  the  electrons 
reside  there.  Therefore,  only  spin-down  electrons  contribute  to  transport  processes.  This 
phenomenon  represents  the  spin-polarization  of  the  system. 

The  distribution  of  a  few  lowest  ground  state  LLs  of  spin-down  electrons,  calculated  from 
the  rate  equation  model  as  a  function  of  magnetic  field,  is  given  in  Fig.  4.  The  decrease  of 
the  LL  population  occurs  at  the  same  magnetic  fields  where  Fermi  energy  drops,  justifying 
the  assumption  that  the  system  is  close  to  equilibrium,  which  was  introduced  to  calculate 
the  concentration  of  spin-up  and  spin-down  electrons. 
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C.  Dark  current 


The  population  in  the  continuum  is  mainly  determined  by  the  transitions  between  the 
ground  LLs  and  the  continuum,  and  the  magnetic  field  influence  on  the  scattering  between 
continuum  LLs  is  of  secondary  importance.  Hence,  the  dark  current  is  determined  by  scat¬ 
tering  from  the  populated  ground  state  LLs,  from  which  the  electrons  are  injected  into  the 
continuum.  Therefore,  scattering  from  the  ground  state  LLs  is  considered  first. 

The  average  scattering  rate  from  the  ground  state  LLs  to  the  continuum  is  defined  as 

rii  E  Wij 

Wa  =  W1 _ 

4  E  nf 

where  i  =  1 1 ,  ,  — 1/2) ,  jt  >  0,  is  a  ground  state  LL,  ns  is  the  doping  density  and  it  is 
assumed  that  the  continuum  LL  population  is  small.  These  are  shown  in  Fig.  5  as  functions 
of  the  magnetic  field,  for  an  electric  field  of  1  kV/cm.  At  a  fixed  electric  field,  an  increasing 
magnetic  field  causes  the  energy  separation  between  any  of  the  ground  state  LLs  except  the 
lowest  ( 1 1,  ji),  ji  >  0),  and  the  continuum  to  decrease  towards  Elo,  significantly  enhancing 
the  corresponding  scattering  rates.  On  the  other  hand,  the  separation  between  the  ground 
state  and  the  continuum  states  increases,  and  becomes  more  distant  from  Elo ,  which  results 
in  a  reduction  of  the  scattering  rate  from  the  j  1, 0)  LL  to  the  continuum.  However,  since 
the  ground  state-continuum  transition  energies  are  remote  from  Elo  in  the  whole  range 
of  magnetic  fields  (a  requirement  in  the  design),  this  decrease  does  not  have  a  large  effect 
on  scattering  from  1 1, 0) .  Since  the  concentration  of  each  ground  state  LL  increases  and 
then  decreases  (apart  from  1 1, 0)  which  increases  and  reaches  saturation),  see  Fig.  4,  all 
average  scattering  rates  increase  until  the  decrease  of  concentration  prevails,  causing  their 
decrease  with  further  increase  of  magnetic  field,  as  shown  in  Fig.  5.  The  described  maxima 
in  1 1 , 1) ,  1 1, 2) ,  1 1,  3)  occur  at  magnetic  fields  of  3.5  T,  1.8  T  and  1.2  T,  respectively,  in 
Fig.  5.  The  other  peaks  in  the  average  scattering  rates  of  1 1,  2)  and  1 1,  3)  in  Fig.  5  appear  as 
the  consequence  of  their  resonances  with  some  of  continuum  LLs,  when  the  scattering  rates 
become  strongly  enhanced  and  prevail  over  the  effects  of  decreasing  population. 

The  dark  current  vs.  magnetic  field  dependence  is  shown  in  Fig.  6.  From  Fig.  5  and  Fig.  6 
one  can  see  that  the  dark  current  follows  the  same  trend  as  the  total  average  scattering  rate 
from  all  the  ground  state  LLs.  Since  only  spin-down  electrons  here  interact  with  phonons, 
the  dark  current  is  spin-down-polarized.  The  zero-field  extrapolation  of  the  results  from  our 
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model  shows  a  reasonable  agreement  with  the  measured  values  of  the  dark  current  in  the  first 
reported  bound-continuum  THz  QWIP31,  as  well  as  in  another,  just  published,  structures32, 
for  the  identical  value  of  electric  held  (1  kV/cm). 

The  dark  current /electric  held  characteristics  at  B  =  3  T  is  shown  in  the  inset  of  Fig.  6. 
Applying  a  stronger  electric  held  causes  an  increase  of  dark  current,  just  as  in  QWIPs 
without  magnetic  held.  This  is  because  the  applied  electric  held  does  not  affect  the  LL 
configuration,  so  despite  the  discretization  of  the  spectrum  this  system  behaves  similarly  to 
the  one  having  simple  parabolic  dispersion  in  the  plane  of  the  layers. 

D.  Responsivity 

The  population  of  higher  index  ground  state  LLs  decreases  with  magnetic  held,  and  hence, 
so  does  the  number  of  photoelectrons  generated  from  them  into  higher  index  continuum 
states  LLs,  see  Fig.  4.  In  contrast,  the  population  of  low  index  ground  state  LLs  and 
of  photoelectrons  generated  from  them  into  low  index  continuum  states  LLs  increase.  The 
contribution  of  higher  index  LLs  photoelectrons  to  the  photocurrent  is  larger  than  that  of  low 
index  LLs  photoelectrons,  because  of  their  larger  transition  rates  to  the  neighboring  periods. 
Therefore,  the  responsivity  generally  exhibits  a  decreasing  dependence  on  the  magnetic  held, 
shown  in  Fig.  7.  However,  this  trend  is  changed  when  an  increasing  number  of  photoclectrons 
in  low  index  continuum  LLs  becomes  more  important  than  the  decreased  number  of  higher 
index  continuum  LLs  photoclectrons,  and  then  the  responsivity  increases.  The  scattering 
rates  from  the  continuum  to  the  ground  state  LLs  increase  with  the  magnetic  held  because 
the  energy  spacing  between  relevant  LLs  approaches  Elo ■  The  responsivity  vs.  magnetic 
held  dependence  is  again  reversed,  reaching  a  maximum  in  Fig.  7,  when  scattering  rates 
from  low  index  continuum  LLs  increase  sufficiently  to  reduce  photoclectron  concentration 
on  them.  The  responsivity  then  continues  to  follow  the  general  trend.  Photocurrent  is  spin- 
down  polarized  because  only  the  spin-down  electron  transitions  occur  via  the  interaction  with 
phonons  and  incident  radiation.  In  comparison  to  previously  reported  bound-continuum 
THz  QWIPs31,32,  at  1  kV/cm  the  order  of  magnitude  of  the  responsivity  in  the  whole  range 
of  magnetic  helds  is  the  same  (~  0.1  A/W). 

The  inset  of  Fig.  7  shows  the  responsivity  as  a  function  of  bias  for  B  =  3  T.  At  a  hxed 
magnetic  held,  the  bias  dependence  of  responsivity  is  qualitatively  the  same  as  in  QWIPs 
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with  no  magnetic  field. 


IV.  CONCLUSION 

A  quantum  mechanical  model  of  electron  transport  in  QWIPs  in  magnetic  field  was 
developed  and  used  to  design  and  analyze  the  performance  of  magnetic  field  tunable  DMS 
THz  QWIPs.  The  tunability  is  based  on  the  magnetically  induced  variation  of  the  band 
offset  potential  profile.  The  QWIP  parameters  were  engineered  to  get  maximal  tunability 
in  the  desired  range  of  the  far  infrared  spectrum,  while  maintaining  low  dark  current  and 
large  responsivity.  The  values  of  current  density  and  responsivity  were  calculated  and  found 
to  be  of  the  same  order  of  magnitude  as  in  GaAs/AlGaAs  THz  QWIPs. 
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Figure  captions 

FIG.  1:  Left-hand  side:  The  conduction  band  profile  of  the  QWIP,  the  energies  of  states 
that  belong  to  one  period  (dashed  lines)  and  the  squared  wavcfunctions  of  the  ground  state 
and  the  continuum  state  with  maximal  oscillator  strength  towards  the  ground  state,  for  the 
electric  field  E  —  1  kV/cm  and  magnetic  held  B  —  5  T.  Right-hand  side:  The  conduction 
band  profile  in  a  magnetic  well  and  in  the  vicinity  of  interfaces:  for  both  spins  at  B  =  0  T 
(solid  line),  for  spin- up  electrons  at  B  —  5  T  (dashed  line)  and  for  spin-down  electrons  at 
B  —  5  T  (dash-dotted  line). 

FIG.  2:  Transition  energies  corresponding  to  maximal  oscillator  strength  for  spin-down 
and  spin- up  electrons  vs.  magnetic  held  dependence,  at  a  hxed  electric  held  of  1  kV/cm. 

FIG.  3:  The  fan-out  of  Landau  levels  from  the  ground  and  the  hrst  continuum  states  for 
both  the  spin-up  and  spin-down  electron  system,  together  with  the  quasi-Fermi  level,  for  a 
hxed  electric  held  of  1  kV / cm. 

FIG.  4:  The  electron  distribution  over  the  hrst  4  ground  state  Landau  levels  vs.  magnetic 
held  dependence  at  a  hxed  electric  held  of  1  kV/cm.  Inset:  the  electron  distributions  over 
the  spin-up  and  spin-down  states  vs.  magnetic  held  dependence  for  the  same  electric  held. 

FIG.  5:  Average  scattering  rates  (see  text  for  explanation)  from  the  hrst  4  ground  state 
Landau  levels  to  the  continuum,  and  the  total  average  scattering  rate  from  all  ground  state 
Landau  levels  vs.  magnetic  held  dependence  at  a  hxed  electric  held  of  1  kV/cm. 

FIG.  6:  Dark  current  density  vs.  magnetic  held  dependence  for  a  hxed  electric  held  of  1 
kV/cm.  Inset:  dark  current  density  vs.  electric  held  dependence  for  a  hxed  magnetic  held 
of  3  T. 

FIG.  7:  Responsivity  vs.  magnetic  held  dependence  for  a  hxed  electric  held  of  1  kV/cm. 
Inset:  responsivity  vs.  electric  held  dependence  for  a  hxed  magnetic  held  of  3  T. 
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